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Abstract 


Probabilistic  depreciation  is  a  method  of  determining  the  proper  depreciation, 
charge  in  each  year  of  an  asset's  service  life,  when  the  service  life  is  a  random 
variable  with  known  distribution.  In  this  paper,  we  discuss  how  the  service 
life  distribution  is  modified  as  we  gain  more  information  about  the  actual  life¬ 
time  of  the  asset.  The  problem  of  determining  the  proper  amount  to  be  charged 
each  year  to  depreciation  while  at  the  same  time  maintaining  the  proper  balance 
in  the  accumulated  depreciation  account  is  considered.  The  analysis  is  done  both 
for  a  single  asset  case  and  for  group  depreciation.  A  final  section  discusses  the 
use  of  Bayesian  analysis  for  estimating  the  particular  form  of  the  service  life 
distribution  while  the  assets  are  in  service. 


1.  Introduction: 

Conventional  depreciation  methods  treat  the  service  life  of  an 
asset  as  a  given  constant.^  When  there  is  some  uncertainty  as  to  the 
actual  service  life,  an  average  (expected)  value  is  used.  The 
depreciation  rates  are  then  calculated  based  on  this  given  service  life. 
Contrary  to  this  deterministic  approach,  a  probabilistic  approach  to 
depreciation  has  been  proposed  in  [2]  where  the  service  life  is  treated 
as  a  random  variable  which  has  a  given  probability  distribution.  In 
probabilistic  depreciation,  the  depreciation  rates  are  calculated  for 
each  possible  service  life  of  the  asset  and  then  the  weighted  average 
is  computed  using  the  service  life  probabilities  as  weights.  Thus  if 
the  asset's  life  is  equally  likely  to  be  1,  2,  or  3  years,  the  depre¬ 
ciation  rate  for  the  first  year  under  the  probabilistic  straight  line 
method  is  the  average  of  100%,  50%,  and  33  1/3  %  or  61  1/9  %.  This  is 
in  contrast  to  the  conventional  straight  line  method  which  computes  an 
average  service  life  first,  2  years  in  this  case,  and  then  calculates 
the  depreciation  rate  as  50%  for  the  first  and  second  year. 

It  was  demonstrated  in  [2]  that  for  single  asset  depreciation, 
conventional  depreciation  methods  result  in  underdepreciation  in  the 
tarller  years  of  an  asset's  service  life.  In  the  case  of  group  depre¬ 
ciation,  conventional  methods  will  typically  result  in  underdepreciation 
throughout  the  service  life  of  the  items  In  the  group.  The  analysis  In 
[2]  computed  the  depreciation  rate  for  each  year  based  only  on  e  single 
estimate  of  the  probability  distribution  of  the  service  life;  namely  the 
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one  known  when  the  asset  was  first  put  into  service.  As  the  asset' 
service  life  expires,  we  obtain  more  information  about  the  actual 
probability  distribution  governing  the  particular  asset  under  study 
and  it  becomes  reasonable  to  modify  the  probability  distribution  to 
include  the  increased  knowledge  that  we  have. 

In  the  numerical  example  above,  if  the  asset  survives  the  first 
year,  we  know  that  the  1  year  service  life  situation  did  not  occur 
and  thus  need  to  consider  only  the  cases  of  2  and  3  years  of  the 
service  life.  Therefore,  the  conditional  probability  distribution  of 
the  service  life,  given  that  the  asset  has  survived  the  first  year, 
is  0.5  for  the  2  year  life  and  0.5  for  the  3  year  life.  Hence, 
the  first  year  depreciation  rate  for  the  asset  after  it  has  survived 
the  first  year  in  service  is  (1/2)  (1/2  +  1/3)  «*  5/12  *  0.42  under  the 
straight  line  method. 

In  this  paper,  we  analyze  how  the  depreciation  patterns  might  be 
affected  by  such  a  "sequential  probabilistic  depreciation"  approach 
as  compared  with  a  "static  probabilistic"  approach  discussed  in  [2J. 
Section  2  discusses  this  problem  for  a  single  asset  case  and  Section  3 
for  group  depredation.  In  Section  4,  an  application  of  Bayesian 
analysis  for  estimating  the  service  life  distribution  is  considered. 

2.  Sequential  Probabilistic  Depreciation;  Single  Asset  Case 

Before  we  proceed  in  our  analysis,  we  shall  define  our  criteria 
in  selecting  depreciation  methods.  In  this  paper,  ve  want  to  set  aside 
the  issue  of  whether  or  not  depreciation  rates  should  be  based  on  the 
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consumption  of  service  potentials,  or  on  the  decline  of  the  market 
values,  or  on  any  other  factors.  We  assume  that  these  economic  consid¬ 
erations  have  been  included  in  the  depreciation  vector  =  (h^,  ^j*  » 

hjj)  where  is  the  proportion  of  the  depreciable  cost  (acquisition 
cost  less  estimated  salvage  value)  to  be  depreciated  in  the  i  -th  year 
if  the  service  life  of  the  asset  is  j  years.  We  shall  assume  that 

>  0  for  all  i  *=  l,2,..,j  and  j  *=  1,2...,  and  that  ^  **1^ j  “  *°  ***** 

the  asset  is  fully  depreciated  by  the  end  of  its  service  life.  For  the 
straight  line  method,  h^  1/j  ;  for  the  sum-of-years-digits  method, 
h^j  *  2(j-i+l)/j(J+l);  and  for  the  double -dec lining  balance  method,  h^  ** 

(1  -  2/j)i"l(2/j). 

With  these  assumptions,  vc  consider  the  situation  in  which  the 
service  life  j  is  not  given  with  certainty  but  is  a  random  variable 
with  a  given  probability  distribution,  the  criterion  used  in  selecting 
a  depreciation  vector  is  that  the  depreciation  rate  under  the  selected 
depreciation  vector  be,  year  for  year,  equal  to  the  expected  value  of 
the  depreciation  rate  under  the  given  probability  distribution.  We  refer 
to  this  property  as  the  unbiasedness  of  the  depredation  rates.  The  method 
of  deriving  unbiased  probabilistic  depreciation  vectors  when  the  asset  Is 
first  placed  in  service  is  discussed  in  [2] • 

In  this  paper,  we  shall  apply  the  sane  criterion  of  unbiasadnes* . 
However,  by  using  the  additional  information  that  is  available  after 
some  number  of  years,  we  shall  try  to  minimise  the  discrepancies  between 
the  ex  ante  depreciation  rates  based  on  the  prior  service  life  distribution 
end  the  ex  post  depreciation  rates  based  on  the  actual  experience  with  the 
asset. 
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Consider  an  asset  with  a  service  life  distribution  given  by  (p^^, 
...,)  where  p^  is  the  probability  that  the  asset  is  retired  at  the  end 
of  the  j  -th  year  in  service.  We  assume  that  retirement  occurs  only 

at  the  end  of  a  year. 

In  order  to  compute  the  unbiased  estimate  of  the  depreciation  rate 
for  the  i  -th  year,  we  must  consider  three  cases.  First,  if  the  asset 
is  retired  in  the  i-1  st  year  or  earlier,  we  have  dj^  =  0  where  d“  is 
the  depreciation  rate  for  the  i  -th  year  when  the  asset  has  been  retired 
in  an  earlier  year.  If  the  asset  is  retired  in  the  i  -th  year,  then 
d°  *  where  d°  is  the  depreciation  rate  for  the  i  -th  year  when  the 
asset  has  been  retired  at  the  end  of  the  i  -th  year.  Note  that  in 
both  of  these  cases,  the  service  life  was  known  with  certainty  so  that 
there  is  no  averaging  required  to  compute  the  proper  depreciation  rate. 

The  third  case  occurs  if  the  asset  survives  through  the  i  -th  year 
so  that  the  service  life  is  known  to  be  greater  than  i.  In  this  situation 
the  probability  that  the  service  life  is  j  (>i)  years  is,  by  the  definition 
of  conditional  probability,  p^/  J  p^  *  Pj/S^  J  ■  1+1,  i+2,.., 

where  ■»  i+1  **  ****  *  Pr*ori  Probability  that  the  asset  is  not 

retired  in  the  first  1  years. 

If,  in  fact,  the  asset  is  retired  in  the  j  -th  year,  the  depreciation 
rate. in  the  i  -th  year  should  be  h^.  Since  ve  do  not,  at  this  stage, 
know  what  the  actual  service  life  will  be,  ve  compute  the  expected  depre- 
eietion  rete  by  everaging  over  all  possible  service  llvos.  '.'hua , 

di  "I J-4+1  hij  <PJ/Si>  "  <l/Si*  Xj“i+lhijPJ 
where  is  the  depreciation  rate  for  the  1  -th  year  when  the  assot  has 

survived  through  the  1  -th  year. 
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Since  whether  or  not  the  asset  is  retired  before,  at  the  end  of, 
or  after  the  i  -th  year  is  not  know,  at  the  tine  the  asset  is  first 
put  into  service,  the  actual  depreciation  rate  for  the  i  -th  year  is 
a  random  variable.  The  expected  value  of  the  actual  i  -th  year 
depreciation  rate  is  given  by 
Vi-l  j"  .  ,o  ,  V  * 

i L  j-l  pj  di  +  pi  di  +  L  j=i+l  pj  di 
=  pihii+  H  j=i+l  Pj  hi  j 
“Ij=iPjhij 

which  is  equal  to  the  static  probabilistic  depreciation  rate  for  the 
i  -th  year  as  derived  in  [2].  Note  that  the  sum  of  these  expected 
depreciation  rates  for  1  *  1,2,..,  la  unity  since 

I  i*l  I  j Pj  hij  "  I  J~1  X  i“l  Pj  hi j  "  I  j»l  Pj  I  i«l  hiJ  ”  1  * 

Suppose  that  the  depreciation  rate  for  the  i  -th  year  when  the 

£  ^ 

asset  has  survived  the  1  -th  year  was  chosen  to  be  d^  )  ,  Then 

the  mean  squared  deviation  between  the  actual  depreciation  rate  and  the 
proper  depreciation  rate  is  given  by 


<l/V  I  j-i+i  pj  ^ij  *  di  >3 


Then, 


<i/si>  l  A+i  ri  Ihu  *  "i*1’  -  <»'v  l 

-o'vtXA+i  p4  Is-m 

-<*/Si)tXA«PjV4i+>l  +  «»+*0*!i5> 

n  K  4  4 

•ince  2,  jni+i  pJ^iJ  “  di^  *  0  by  tlt*  chotce  • 
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Thus,  the  mean  squared  deviation  b  tween  the  actual  depreciation 

rate  and  the  proper  depreciation  rate  is  minimized  when  we  choose 

*  + 

4i  ■  «l  ■ 

So  far  we  have  concentrated  on  computing  the  actual  depreciation 
rate  for  the  i  -th  year  (i  =  1,2,..,)  to  give  as  accurate  a  figure  as 
possible  for  use  in  the  income  statement.  We  now  consider  what  the 
accumulated  depreciation  rate  should  be  at  the  end  of  the  1  -th  year 
after  the  depreciation  rate  for  the  i  -th  year  has  been  charged. 

Let  a^  be  the  proportion  of  the  depreciable  cost  that  has  already 
been  depreciated  by  the  end  of  the  1  -th  year.  Clearly,  if  the  asset 
has  been  retired  by  the  end  of  the  i  -th  year,  a^  should  equal  1  . 

If  the  asset  is  still  in  service  at  the  end  of  the  1  -th  year,  a^ 
should  equal  1  minus  the  sum  of  the  expected  depreciation  rates  in 
each  of  the  remaining  years  in  order  for  a^  to  be  an  unbiased  estimate 
of  the  proper  accumulated  depredation  rate  when  the  asset  is  still  in 
service  at  the  end  of  the  i  -th  year. 

The  expected  value  of  the  depreciation  rate  in  the  k  -th  year  given 
that  the  asset  is  in  service  at  the  end  of  the  i  -th  year  is 

(‘'V'Pj-M  "#  +  ’V  *Vk  +  E  J*k+1  PJ  * 

‘  <‘*1> X£  \,  ' 

Therefore,  the  sum  of  the  expected  depreciation  rates  in  all  years  after 
the  i  -th  is 

*l/V  Zk-l+llj-k  PJ  *Vj  Zfc^I+i  *kj  * 

*«*•.  "  1  *  U/St>  £  JHL41  ZiU+1 

"  <l/V  Xj-4+i  , 

•iace  *7^*1  “l  * 
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Note  that  a*  equals  the  sun  of  the  expected  depreciation  rates  that  should 
have  been  taken  in  years  1  to  1  if  it  were  known  that  the  asset  would 
still  be  in  service  at  the  end  of  the  i  -th  year. 

Knowing  a*  ,  we  can  compute  the  expected  accumulated  depreciation 
rate  after  i  years  to  be 

L  j«l  Pj  *  1  +  (I  J*4+l  Pj5  ®i 

“I  j=l  PJ  +  Si(1/Si)  I  j=i+l  pj  Xk-1  \j 

~  I  j~l  Pj  /!  k=l  hkj  +  I  j*i+l  Pj  I  k=l  \i 

*Xk~l  I!  j*k  Pj  \j  +  T,k®l  T  j=ivl  pj  Ntj 

“  I  k«l  Ij*%  Pj  \j 

which  is  exactly  equal  to  the  accumulated  depreciation  rate  under  the 
static  probabilistic  method  as  derived  in  [2).  Since  this  quantity, 
which  is  less  than  unity,  is  a  convex  combination  of  1  and  a^  ,  we 
know  that  the  accumulated  depreciation  for  nonretired  assets  Is  always 
less  using  sequential  probabilistic  depreciation  (a*)  than  when  using 
static  probabilistic  depreciation. 

Since  the  accumulated  depreciation  rate  and  depreciation  rates 
for  each  year  were  derived  separately,  it  is  not,  in  general,  true 
that  the  accumulated  depreciation  rate  is  the  sum  of  depreciation  rates 
in  each  year.  If  the  asset  is  still  in  service  at  the  end  of  year  1*1, 
Its  accumulated  depreciation  rate  is  given  by  a^j  .  f  it  survives 
the  i  -th  year,  the  depreciation  rate  is  given  by  d*  .  The  difference, 

+  d^)  ,  may  be  considered  as  an  entry  to  an  account  for 
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adjusting  prior  income,  and  is  attributable  to  the  uncertainty  in  the 

service  life  of  the  asset.  Similarly,  if  the  asset  is  retired  in  the 

i  -th  year  ,  the  depreciation  rate  is  d°  and  the  adjustment  is 

given  by  1  -  (a*^  + 

An  interesting  feature  of  this  method  is  that  at  does  not  neces¬ 
sarily  increase  as  i  increases*  For  example,  if  p^  —  p^q  =  0.5 
(all  other  =  0)  ,  and  if  straight  line  depreciation  is  used,  the 
accumulated  depreciation  rate  after  one  year,  a*  ,  is  .5(1/2)  +  .5(1/10) 

=  3/10  .  If  the  asset  is  still  in  service  at  the  end  of  the  second  year, 
its  accumulated  depreciation  rate  a^  shou!.d  be  1*(2/10)  ,  a  decrease  of 
1/10  over  a*  .  Since  the  asset  is  still  in  service  at  the  end  of  the 

second  year,  i.e. ,  the  asset's  service  life  is  now  known  to  be  10  years 

with  certainty,  the  first  year's  depreciation  rate  should  have  been 
1/10.  This  plus  the  depreciation  rate  in  the  second  year  of  1/10  equals 

X  X 

*2*  The  difference  between  a^^  =  3/10  and  what  it  should  have  been, 

1/10  ,  is  attributable  to  the  uncertainty  in  the  service  life.  Thus  a 
proper  accounting  treatment  is  to  depreciate  1/10  in  the  second  year  and 
stake  a  debit  entry  of  2/10  to  the  accumulated  depreciation  rate,  crediting 
to  an  account  for  adjusting  prior  income. 

In  general |  the  expected  value  of  the  adjustment  required  in  the 
1  -th  year  Is 

U*  -  (•*.,  +  <£>]  St  +  U  -  +  d°)l  Pt 

"Ij-i+1  Pj  Z.W  ^j  *  Ij*l  hij  Pj  *  Pi  “  Ij"*  Pj  Ik=l 

"I  jK+l  Pj  Ik»l  ^j  *  I  j-i  Pj  Ik=l  Ncj  +  Pi 

-  0 
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so  that  the  adjustment  is  due  solely  to  the  inherent  fluctuations  caused 
by  the  variability  in  the  service  life  of  the  asset. 

3.  Sequential  Group  Depreciation 

The  sequential  probabilistic  depreciation  method  developed  in 
Section  2  can  be  easily  applied  to  group  depreciation  in  which  assets 
of  a  similar  nature  are  grouped  together  and  depreciated  in  a  single 
account.  At  the  end  of  any  year,  the  original  group  of  assets  can  be 
partitioned  into  three  groups;  those  that  were  retired  in  year  i-1  or 
earlier ,  those  that  have  just  been  retired  in  year  i  ,  and  those  that 
are  still  in  service.  The  items  in  the  first  group  have  already  been 
fully  depreciated  and  no  further  charges  are  required.  Those  items 
which  have  just  been  retired  are  charged  at  a  rate  equal  to  d°  =  h^ 
while  those  items  still  itf  service  are  charged  at  a  rate  equal  to 

di  *  (1/Si>  Lj=i+1  hij  Pj  ' 

Let  N  be  the  number  of  items  originally  placed  in  service  at  the 
start  of  the  first  year.  Let  be  the  number  of  items  that  have  been 
retired  at  the  end  of  the  i  -th  year  and  let  be  the  number  of  items 
that  were  retired  just  at  the  end  of  the  i  -th  year.  Thus, 

ni  “Vl  +  wi  * 

Therefore  the  i  -th  year  depreciation  rate,  based  on  the  total 
depreciable  cost  for  the  N  items,  is  given  by  » 

ni-l  dI  +  d°  +  (N  -  n^)  d* 

TH  N~ 

“  hii  +  (N  '  V  r  « 

r  4  nst-1  Ij^ihijPj  • 


where  d^  ,  d°  ,  and  d*  are  as  defined  in  Section  2.  In  order  to  compute 
the  expected  depreciation  rate  in  the  i  -th  year,  we  need  to  compute 
only  E(m^/N)  and  E(N  -  n^/N)  since  d°  and  d*  are  fixed  positive  numbers 
and  =  0  .  The  probability  mass  functions  for  the  random  variables 
and  N  -  n^  are  given  by  a  binomial  distribution  with  parameters 
n  (for  both  variables)  and  p  =  p^  and  respectively.  Since  the 
expected  value  of  a  binomial  distribution  with  parameters  n  and  p  is 
np  we  have  that 

EOn^/N)  ~  p^  and 

E(N  -  n^)/N  *  S± 

Therefore  the  expected  i  -th  year  sequential  group  depreciation  charge  is 
given  by 

tpt  t>uJ  +  »t  (1/Sj)  J  j=i+i  hij  pj 

the  static  group  probabilistic  rate  computed  in  [2J. 

A  similar  analysis  shows  that  the  accumulated  depreciation  rate 
after  1  years  is 

"t  • 1  +  (tl  •  at  ■=nt  I  ji+i  pj  Iw  \j  • 

T  s  T  1s^~ 

Therefore,  the  expected  accumulated  depreciation  rate  after  i  years  is 

“  I  j»l  +  I  j»i+l  Pj  Ik»l  hkj 

the  sccurnulated  depreciation  rate  after  i  years  under  the  static  probabilistic, 
method  (see  [2)). 
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As  In  the  single  asset  case,  an  adjustment  is  generally  needed  each 
year  to  keep  the  accounts  balanced.  The  amount  of  this  adjustment  is 

(1  -  hu>  [m.  -  Pi(H  -  Vl)l  *1^  Pj  ££  ly  &  -  IVl] 

NS^  Kbi  NSi-l 

and  one  can  readily  check  that  the  expected  value  of  this  adjustment  is 
zero. 


4.  Bayesian  Analysis 

For  the  calculations  described  in  the  previous  sections,  we  have 
implicitly  assumed  that  the  service  life  distribution, {p  ;  j  =  1,2,..,], 
was  precisely  known.  In  [2]  we  discussed  the  use  of  parametric  distri¬ 
butions  such  as  the  normal,  rectangular,  Poisson,  and  geometric  to  provide 
convenient  representations  for  the  service  life  distributions.  In  using 
these  distributions,  one  makes  the  assumption  that  the  probabilities 
of  asset  retirement  follow  about  the  pattern  indicated  by  the  distri¬ 
bution  so  that  only  one  or  two  parameters  need  to  be  specified  in  order 
to  completely  characterize  the  distribution. 

In  practice,  one  may  be  willing  to  accept  the  assumption  that  a 
particular  distribution  is  appropriate  and  yet  still  be  unsure  about 

the  parameters  that  characterize  the  distribution.  In  this  case, 

3/ 

Bayesian  analysis  might  prove  useful.- 

Assume  we  have  a  group  of  N  identical  assets  each  having  an  independent 

and  identically  distributed  service  life  distribution  given  by  [pj  ;  j«=l , 2 , . , , } 

Then,  the  probability  that  exactly  n^  are  retired  in  the  first  year  (assuming 

p^  >  0)  is  given  by 

,Nv  n,  vN-n. 

<n)  pi  i  a  -  ?i>  i  • 
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Assume  also  that  is  not  known  for  certain.  The  question  is  how  we 
can  use  our  prior  information  on  p^  with  the  actual  retirement  number, 
say  ,  in  the  first  year  to  get  &  new  estimate  on  and  indirectly 
on  P2>P3»**>  •  We  could  then  continue  in  this  manner  at  the  end  of 
each  year  to  note  the  retirement  pattern  and  update  our  estimate  of  the 
parameter (s)  of  the  distribution  under  study. 

Let  f^(p^)  be  our  prior  density  function  on  p^  where 
f^(p^)  >  0  for  0  <  <  1  and 

Jo  W  dP  “ 1  • 

After  observing  n^  retirements  in  the  first  year,  our  new  (posterior) 

density  function  on  p^  (denoted  by  f^(P^)  )  is  given  by 
n.  N-n, 


W  " 


'(l-Pj) 


Lf(Pl) 


fl  n.  N-n 

Jo  -i 


0  <  p  <  1 


f(pL)  dp 


In  general,  however,  one  would  not  be  solely  interested  in  but  rather 
in  what  the  information  on  p^  implies  about  the  parameter (s)  that  are 
specifying  the  complete  service  life  distribution. 

For  example,  suppose  one  feels  that  the  geometric  distribution  is 
•  reasonable  approximation  to  the  service  life  distribution.  Therefore 

»j  ‘  (1  -  •  3  "  l'2 .  y 

One  way  of  estimating  p  is  by  noting  that  the  expected  lifetime  is  given 

b)r  1  )’i 3  pj  * l/p  .  However,  one  may  not  have  had  sufficient 

experience  with  items  sirjiliar  to  the  assets  currently  being  used  to 

place  complete  faith  in  a  point  estimate  of  p  .  Therefore  one  might  wish 

to  select  a  prior  density  on  p  and  let  the  experience  with  the  ectual  items 

now  in  service  also  influence  our  estimate  of  p  . 


For  exposition  purposes,  it  i£  convenient  in  this  situation  to 
choose  a  beta  cistribution  as  the  prior  density  for  ;  i .e. , 

f(p)  *  A0  Pr“1(l  -  pf'1?  »  0  <  p  <  1;0  <  r  <  t  , 


where  = 


(t-l)i _ 

(r-l) !  (t-r-1) ! 


The  parameters  r  and  t  may  be  estimated  by  using  the  relations  that 
E(p)  =  r/t  and 
Var(p)  =  r/t2  . 

Another  interpretation  of  how  r  and  t  may  be  estimated  will  become  clearer 
as  we  proceed. 

•  Therefore,  if  ^  out  of  N  items  are  retired  in  the  first  year  in  service 


fx(p) 


"l  N'ni  r-l  t‘r“1 

A0  P  l(l-p)  p  A(l-p) 

fn,  N-n  r-l  t-r-1 
AnJp  i(l-p)  1  p  (1-p)  dp 


(N+t-1)! 

(r+nj-1)!  (N+t-r-n^-1)! 


r4n.-l  N+t-(r+n  )-l 
P  1  <1-P>  1 


which  is  a  0-  distribution  with  parameters  r-Hij,  ,  N+t  . 

At  the  start  of  the  second  period,  we  have  N  -  n^  items  still  in  service 

Suppose  that  nj  -  n^  (>  0)  items  fail  in  the  second  period.  The  probability 

of  this  event,  given  p  ,  is 

„  n.-n.  N-n, 

(N-ni  }  p  2  i  (l.p)  2  # 

V"l 

Therefore  if  our  prior  on  p  is  given  by  f^(p)  above,  the  posterior  (after 
the  second  period)  density  of  p  is 
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iu-n,  N-n-  r+n.-l  N+t-Cr+n^-l 
tAP)=koP*  \l-P)  p  (1-p)  1 


=  a2  p 


r+c2~l 


N+(N-n1)+t-(n24x)-l 


where 


A2  - 


(1-P) 

(N+N-n^+t-1)! 


(r+n2-l)I  (2N+t-n1-r-n2-l)l 
Proceeding  in  this  manner,  it  can  be  readily  shovm  by  induction  that  if 
-  nj-l  *tems  have  been  retired  in  the  j-th  period  (j  —  l,2,..,k;  n^  —  0)  , 
the  posterior  density  of  p  after  the  k-th  period  is  given  by 


.  <»-!>  '  »' _ 


where 


(r-K^-1)!  (t+££^  (N-np  -  (r+i^)  -  1)! 


Mote  that  f^(p)  d pends  on  the  entire  history  of  retirements  (as  represented 
by  n^, ..,n^)  and  not  just  on  ,  the  cumulative  number  of  retirements  at 
the  end  of  the  k-th  year. 

One  can  see  that  the  initial  parameter*  r  and  t  can  be  interpreted  as 
if  one  had  a  prior  experience  of  observing  t  of  these  identical  item9  in  use 
«t  the  start  of  a  period  with  r  of  them  having  been  retired  by  the  end  of 
the  period.  Thus  the  greater  t  is,  the  larger  the  influence  of  the  prior 
density  on  the  posterior  densities  of  later  periods.  In  the  limit,  where 
we  know  p  with  certainty  so  that  a  Bayesian  analysis  is  not  useful,  r  and  t 
both  go  to  infinity  with  the  ratio  r/t  remaining  equal  to  p  • 

Finally,  to  compute  the  probabilistic  depreciation  rate  and  accumulated 
depreciation  rate  after  the  k-th  period,  we  need  eo  determine  Pj  for  j  » 
k+l,k+2,..,  .  We  know  that  for  a  given  value  of  p  ,  p^jp  *  (i-p)  ^  ^  ^p  . 
therefore,  to  find  the  unconditional  value  of  p^  ,  we  multiply  by  the  density 


function  for  p  ,  f^(p)  aR^  integrate  over  0  to  1  ;  i.e. , 


Pj  -  J  J  (PjIp)  fk(p)  dP  =  aJ  J  pt*k<i-P)r‘+Ii4<,i-"i)-<rrtk)+«-k-l>-1'dp 

ItfgS  (N-n^-Ut 


<r4V 


[t+Zl4  (N-n^-Oytx)-!]! 


for  j  =  k+l,k+2, 


Thus,  for  example, 
r4pk 

Pk+l  •  r  k-i  . 
^1=0  (N“ni) 


t+L  t*o  (N-ni)-(nk+T> 


Pk+2  ‘+ii4  <K-ni)  £+Im  <“-“i>+i 

•o  that  if  k  «*  l  (have  observed  one  period  with  failures),  our  revised 
estimates  of  p2  end  p^  are  given  by 

p  .nd  D-  *+°l 

t-W  t+K  t-HI+1 

(One  can  easily  verify  that  the  prior  estimates  of  p2  and  p^  are  given  by 

p,  ■  r  and  p«  -  r  .) 

*  t  1  tt+1 


The  values  of  p^  for  j  “  k+1  ,k+2 , , . ,  are  then  Inserted  into  the  formulas 
developed  in  Section  3  to  compute  the  k-th  year  depreciation  rate  and 
accumulated  depreciation  rate  under  the  sequential  Probabilistic  method. 

In  the  above  example,  the  posterior  was  always  a  %-  distribution  with  the 
two  parameters  a  simple  linear  f<  etiou  of  the  prior  estimates  r  and  t  , 
and  the  observed  retirement  history.  The  ease  with  which  we  obtained  the 


posterior  density  of  p  given  n^,...,nk  was  due  to  the  special  structure  that 
was  assumed:  the  gcmetrfc  distribution  of  the  sat  vice  lives  and  the  J-  distrl 
button  as  the  prior  density  on  p  .  Departures  from  either  one  of  these  asaump- 
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tions  will  result  in  considerably  more  complication  in  the  computation  of 
the  posterior  density  function. 

For  example,  in  [2]  we  postulated  the  use  of  a  Poisson  distribution 
as  a  good  approximation  to  actual  mortality  curves  of  equipment.  With  the 
Poisson  distribution,  we  define  p^  -  [X*  V(i-l)!]  e  ^  for  i  -  1,2,..,  . 

We  would  like  to  be  able  to  treat  the  case  in  which  X  is  a  random  variable 
specified  initially  by  a  prior  density  function  f(X)  .  The  authors  are 
unaware  of  any  prior  density  for  X  that  after  observing  a  sequence  of 
retirements  will  yield  a  posterior  density  for  X  of  the  same  general  form. 
However,  if  a  prior  density  on  X  is  chosen,  it  is  still  possible  to  carry 
out  the  required  computations  to  determine,  numerically,  the  posterior 
density  of  X  . 

Let  X  be  specified  by  a  prior  density,  f(X)  .  Then  with  N  identical 

items,  the  probability  that  are  retired  in  the  first  period  is  given  by 

N  .  „  .  N-n. 

|H,X]  «  (  i^Xe'Vl  (l  -  e"*)  1  . 

Therefore  after  observing  n^  retirements,  the  posterior  density  of  X  1* 
ict 

fj(X)  "A^(e"*)  l(l  -  e“k)  lf(X)  where  is  chosen  so  that  f^(X)dX  *  1  . 

In  general,  let  p^(X)  be  the  probability,  given  X  ,  that  ah  item  is 
retired  in  the  k  -th  period  given  that  it  Is  still  in  service  at  the  end 
of  the  k-1  et  period.  Therefore, 

^(X)-  tX^/fk-Dt]/^*^  XJ’l/(J-DI  J  . 

Also,  let  fj^(X)  be  the  posterior  density  of  X  based  on  the  actual  retire¬ 
ment  experience  up  through  the  k-1  et  period.  Then  if  N  -  iters  are 
still  in  service  et  the  start  of  the  k  -th  period  and  n^  -  retire- 
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oients  are  observed  in  the  k  -th  period,  the  new  posterior  density  of  X 
is  given  by 

fk(x)  -  Aklpk<x)]*‘k  ^  11  "  Pk(x)]  N  ^  fk-l(X)  for  x  >  0  ;  k  55  » 

with  chosen  so  that  J  *  f^O.)  d\  =  1  . 


Thus  to  compute  the  sequential  depreciation  rate  for  the  k  -th  period, 
we  need  to  determine  the  unconditional  value  of  for  j  =  k+l,k+2,..., 
as  given  by 


Jo 


[XJ‘l/(j-l)!J 


Ei=k+1  fX  /(f-DH 


fk(X)  dX  for  j  =  k+l.k+2. 


The  above  recursive  computations  are  straightforward,  albeit  tedious,  and 
can  be  programmed  for  solution  on  a  digital  computer. 


mmm 


FOOTNOTES 


1.  For  a  comprehensive  treatment  of  conventional  depreciation  methods,  see 
Grant  and  Norton  [1]  . 

2.  If  the  asset  is  retired  in  one  of  the  first  few  years  of  its  possible 
service  life,  this  adjustment  might  be  quite  substantial.  However,  to 
offset  this  effect,  the  salvage  value  may  be  larger  than  expected  because  of 
the  few  years  in  service  and  hence  may  reduce  the  amount  of  this  writeoff. 

3.  See  Raiffa  and  Schlaifer  {3]  . 

4.  A  more  general  form  of  the  geometric  distribution  might  allow  for  a 

minimum  lifetime  of  m  years  so  that  *  0  for  J  *  l,2,..,m  ,  and 
Pj  ■  (1  -  p)^*’m*’1p  for  j  ■  mf  l,mf2, . . . ,  . 

5.  The  beta  distribution  is  the  conjugate  prior  (see  (3]>  for  the  Bernoulli 
process  being  described  here  where  the  probability  of  retirement  in  each  year 
is  •  constant,  p  ,  for  those  items  which  are  still  in  service. 
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x  Probabilistic  depreciation  is  a  method  of  determining  the  proper  deprecia¬ 
tion  charge  in  each  year  of  an  asset's  service  life,  when  the  service  life  is  a 
random  variable  with  known  distribution.  In  this  paper,. we  discuss  how  the 
service  life  distribution  is  modified  as  we  gain  more  information  about  the 
actual  lifetime  of  the  asset.  The  problem  of  determining  the  propel  amount  to 
be  charged  each  year  to  depreciation  while  at  the  same  time  maintaining  the 
proper  balance  in  the  accumulated  depreciation  account  is  considered.  The  analysis 
Is  done  both  for  a  single  asset  case  and  for  group  depreciation.  A  final  section 
discusses  the  use  of  Bayesian  analysis  for  estimating  the  particular  form  c*f 
the  service  life  distribution  while  the  assets  are  in  service. 
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